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Please replace the paragraph beginning at page 23, line 18 with the 
following rewritten paragraph: 


The implicit Euler integration method is illustrated in the flow chart of 
Fig. 6 for the vector function y = f(y,t) (where y = (q,u), q representing the position 
states and u the velocity states of the molecular system). The function / includes both 
the multibody system dynamics and the forces such as electrostatic attraction and 
repulsion, van der Waal's forces, and solvation forces. After an entry step 79, the first 
operation step 80 updates the Iteration matrix G . For all implicit integration methods, 
the Iteration matrix G has the form G = I-aJ , where / is the identity matrix, a is 
some scalar function of the timestep h n , the timestep between time t n and t n _ x , and J , 

the Jacobian given by J = — . For the implicit Euler method, a = h n . In passing, for 

dy 

additional savings in computer time, it should be noted that a very efficient method of 
computing Jacobian matrices from the residual form of equations is covered in previously 
cited co-pending U.S. Patent Appln. No. 10/053,348, entitled "METHOD FOR 
ANALYTICAL JACOBIAN COMPUTATION IN MOLECULAR MODELING," filed 
of even date and is assigned to the present assignee. As in the case of the present 
invention, the same referenced patent application also describes the use of internal 
coordinates to describe the state of the molecular system. For example, the rotation of 
one part of the molecule is described with respect to another part, rather than with respect 
to an external referenced coordinates. This further increases computing efficiency. 

REMARKS 

The second preliminary amendment corrects two typographical errors in 
mathematical symbols made in the last preliminary amendment. For the convenience of 
the Examiner , the chan g es in the VERSION WITH MARKINGS TO SHOW CHANGES 
MADE are highlighted in yellow. 


